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Abstract 

We propose a new framework called Evolving Systems to describe the self-assembly, or 
autonomous assembly, of actively controlled dynamical subsystems into an Evolved System with 
a higher purpose. Autonomous assembly of large, complex flexible structures in space is a target 
application for Evolving Systems. A critical requirement for autonomous assembling structures 
is that they remain stable during and after assembly. The fundamental topic of inheritance of 
stability, dissipativity, and passivity in Evolving Systems is the primary focus of this research. 
In this paper, we develop an adaptive key component controller to restore stability in Nonlinear 
Evolving Systems that would otherwise fail to inherit the stability traits of their components. We 
provide sufficient conditions for the use of this novel control method and demonstrate its use on 
an illustrative example. 

Keywords: Evolving Systems, stability, adaptive control, dissipativity, passivity 


1. Introduction 

We have proposed a new framework called Evolving Systems to describe the self-assembly, 
or autonomous assembly, of actively controlled dynamical subsystems into an Evolved System 
with a higher purpose [1,2], The components of an Evolving System self-assemble, or mate, to 
complete the Evolved System. The Evolving Systems framework provides a scalable, modular 
architecture to model and analyze the subsystem components, their connections to other compo- 
nents, and the Evolved System. Ultimately, once all the components of an Evolving System have 
joined together to form the fully Evolved System, it will have a new, higher purpose that could 
not have been achieved by the individual components collectively. 

Autonomous assembly of large, complex structures in space, or on-orbit assembly, is an 
excellent application area for Evolving Systems. Future space missions will require on-orbit 
assembly of large aperture (greater than 10 meters) space systems, possibly at distant locations 
that prohibit astronaut intervention [3], Other applications for Evolving Systems can be found 
in [2, 4]. These applications motivate the development of flexible structure Evolving Systems, 
which are mechanical dynamical systems consisting of actively controlled flexible structure com- 
ponents joined together by compliant forces. 
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A critical requirement for autonomous self-assembling structures is that they remain stable 
during and after assembly. We say that a subsystem trait is inherited by an Evolving System 
when the system retains the properties of the trait after assembly. The inheritance of subsystem 
traits, such as controllability, observability, stability, passivity, and dissipativity, is an important 
research topic. 

The fundamental topic of inheritance of stability, passivity, and dissipativity in Evolving 
Systems is our primary focus. In this paper, we develop an adaptive key component controller to 
restore stability to Nonlinear Evolving Systems that would otherwise fail to inherit the stability 
traits of their components. We provide sufficient conditions for the use of the control method. 
Sufficient conditions for the inheritance of passivity and dissipativity traits in Nonlinear Evolving 
Systems are also given. For linear time-invariant (LTI) systems, we show that strict dissipativity, 
strict passivity, and strict positive real are equivalent properties. Finally, we demonstrate the 
adaptive key component stability restoring controller on an illustrative example. 

Decentralized control theory and analysis have been applied to the control of large inter- 
connected systems; see the excellent survey paper by Nils Sandell [5] on this topic. Generally, 
decentralized control has been used to decrease the complexity of the control issues affecting 
large interconnected systems. A decentralized adaptive control approach presented in [6] differs 
from the approach developed here in many ways, including its use of the system’s inputs and 
outputs. The ideas addressed by decentralized control theory are related, but not equivalent, to 
the Evolving Systems viewpoint. 

On the experimental side, a research group at the Information Sciences Institute at the Uni- 
versity of Southern California (USC) has been conducting research in self-reconfigurable, au- 
tonomous robots and systems. They have conducted experimental work to study the feasibility 
of techniques for assembling large space structures as part of their FIMER (Free-flying Intel- 
ligent MatchmakER robots) project [7, 8], This group uses a distributed control method with 
simple proportional derivative control laws for the self-assembly of components. 

2. Mathematical Formulation of Evolving Systems 

In this section we give the general mathematical formulation of nonlinear time-invariant 
Evolving Systems. See [2, 4] for a more detailed description of the Evolving Systems formula- 
tion. Often, applications have significant nonlinearities that must be accounted for in the system 
model. Hence, the framework presented here allows nonlinear dynamics in the system model. 

The components are the building blocks of the Evolving System. Consider a system of L 
individually actively controlled components, where the dynamical equations of the components 
are given by 


f = f(\„ u,); x,(0) = x i0 
1 y i = g/(x/,u,) 


where i = 1 , 2, . . . , L, x, = [x'j x' 2 • • • x‘ n . j is the component state vector, 

x, = |xj x' 2 • • • xj, | , u; = \u\ «!, • • • u' m ] is the control input vector, 

y ; s j}’' y 2 • • • y' p _ | is the vector of sensed outputs, and x, 0 = [x' lo x^ • • • xJ, ;o ] is the 

initial condition. Note that n,- is the dimension of the state vector x,-, m,- is the dimension of the 
control vector u,, and /;, is the dimension of the output vector y,. 
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Each component uses local control to remain stable and to meet its performance require- 
ments. Local control means dependence only on local state or local output information, i.e., 
u, = hj(Xi) or u, = hj(yi). In general, the local controller for an Evolving System component 
would have the form given by 


u, = hj(Zi ) 

z.j = lj( Z/, y ,) 


( 2 ) 


where h, and /, are control operators and z represents the dynamical part of the control law. 

A key concept in Evolving Systems is the evolution of the connection between components, 
which is controlled by a connection parameter, e, that multiplies the forces between components. 
It is critical to note that connection parameters evolve independent of time. We are ignoring 
time in our formulation because it is assumed that joining of components is not time critical and 
because we are interested in studying the system as the components mate. In this framework, the 
evolution of the connection parameter comprises the homotopy 0 < e < 1 , not just the endpoints 
where e — 0 or e — 1 . Components are unconnected when e = 0, and connected when e = 1 . An 
Evolving System is fully evolved when all of the connection parameters joining the subsystem 
components equal 1 . 

Compliant forces join the components, so no degrees of freedom are lost as a consequence 
of two components connecting in a rigid manner. In flexible structure Evolving Systems, the 
compliant forces are usually springs. 

The subsystem components of the Evolving System with interconnections included have the 
form 


X,- = f(xi, u,) + ^ Ljk/jix, u) 

7=1 

y; = g<(x„ u,) 


(3) 


where x = [xj xi ••• X/ 1 , u = |ui 112 ••• U/J , and the function kjfx, u) represents 

the interconnection between components i and j with connection parameter 0 < e, ; - < 1 . 

When the subsystem components join to form an Evolved System, the new entity becomes 


X = /(X, u) 

y = t ?(x, u) 


(4) 


The above is the most general form of an Evolving System. Next, we introduce the framework 
for Evolving Systems of flexible structures. 

A large and important class of nonlinear systems are linear in their input. In [2], we developed 
a state space description of a finite element model representation of a flexible structure Nonlinear 
Evolving System that is linear in its inputs. State space equations for individual components of 
an Evolving System, including component connections, are given by 


x, = A,-(Xj) + B,'(x, )u, + e,yA, 7 (x); x,( 0 ) - x, c 

7=1 

y i = C;(Xj) 


where i = 1 , 2, . . . , L, x, = x' 2 ■ ■ ■ x‘ n . j is the component state vector, 

u ' m , J is the control input vector, y ,■ = [y'j 
3 


u, s u\ 


y 2 


’ P/ 1 is the 


(5) 


vector 



of sensed outputs, x, 0 is the initial condition, A ,(x,), B,(x,), and C,(x,) are matrices that are 
functions of the state and have dimension n,- x n„ n, x m„ and /;, x n,-, respectively, A ;/ (x) is the 
connection matrix, and 0 < e, ; < 1 is the connection parameter. The connection matrix, A i; (x), 
has dimension n,- by dim(x), where n, is the dimension of the state vector x,- corresponding to 

L 

component i and dim(x) = ^ The notation (A(x). B(x), C(x)) is used to denote an Evolving 

k= 1 

System represented by (5). 

3. Impedance-Admittance Formulation of Contact Dynamics 

In this section, we formulate the contact dynamics in Evolving Systems in terms of me- 
chanical impedance and admittance, as first described in [9]. For many dynamical systems, the 
impedance-admittance form is a useful tool for modeling the contact dynamics of components 
[10], The impedance of a mechanical system is determined by the equation / = Z(v), where / 
is the force exerted by the system, v is the velocity of the system, and Z is the impedance of the 
system. The admittance of a mechanical system, Y, is determined by the equation v = Y(f) and 
it is the inverse of the impedance of the system, e.g., Y = Z 1 . Impedance and admittance can 
be seen as nonlinear operators describing the relationship between the output of a mechanical 
system, or the force it exerts at a contact point, with the input of the system, or the velocity at the 
contact point. 

When two components join at a point of contact, their velocities are equal and the forces 
exerted are equal and opposite. Thus, we can model the contact dynamics of two components 
mating in an Evolving System as the feedback connection of the impedance of one component 
with the admittance of the other component. We now introduce two nonlinear operators Y\ and 
Z 2 that provide the admittance and impedance formulation of the contact dynamics of Nonlinear 
Evolving Systems components. These operators relate the force and velocity at the contact point 
of two mating components as given by the equations Vi = Y\ (fj ) and f 2 = Z 2 (v 2 ). In linear 
time-invariant systems, these operators can be easily calculated using Laplace transforms. For 
nonlinear components, the admittance and impedance operators cannot be easily found. How- 
ever, this does not invalidate the analysis provided here, which provides a foundation for adaptive 
key component control and inheritance of passivity and dissipativity traits in Evolving Systems. 

We can use the nonlinear impedance and admittance operators to obtain a state space rep- 
resentation of Evolving System components in admittance-impedance form. Following is a two 
component Evolving System in admittance-impedance form 

f X,- = A,(x,) + eB,-(x/)u/ + Bf (x,)uf 

\ y ,• = C,-(x;) (6) 

1 yf = c?(x,) 

where i = 1,2, A,(x,) represents a component with local control, interconnections between com- 
ponents are represented in B,(x,)u,, control augmentation for stability of the joined system is 
accomplished through B A (x ( )u A , and 0 < e < 1. The Evolved System is the feedback connection 
of the two subsystems, i.e., where (yi = U 2 ) and fui = — y 2 ). 

4. Adaptive Key Component Stability Restoring Controllers 

The Evolving Systems framework facilitates flexibility in the design of systems from com- 
ponents. Stability of an Evolving System can be lost when stable components assemble. There 
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are many advantages, such as cost and reliability, to enable assembly of Evolving System com- 
ponents with little or no modification to the individual components. In this section, we introduce 
adaptive key component control to restore stability to an Evolving System by augmenting the 
controller on a single subsystem component, using only the input-output ports on that compo- 
nent. A fixed gain key component controller was first proposed in [11]. 

Aerospace systems are often difficult and costly to model due to their complexity and their 
uncertain operating environments, making them ideal candidates for adaptive control, which does 
not require detailed knowledge of the system parameters. Adaptive key component control was 
first proposed in [12]. In this approach, one component is chosen to be the key component, 
which is then augmented with a direct adaptive control law that only uses inputs and outputs 
available on that component. Since the key component controller is restoring stability to the 
system without using inputs or outputs on any of the other components in the Evolving System, 
the system needs to be controllable and observable from the key component. 

Now we give equations for an Evolving System with a key component controller. Without 
loss of generality, we can let component 1 be the key component. We may think of component 2 
as being the rest of the Evolving System to which the key component and its adaptive controller 
are connected. Consider a two component Nonlinear Evolving System in impedance-admittance 
form given by 


f X,- = A;(x,) + eB,(x,)u ; - + Bf (x,)uf 

j y,- = C,(x,) (7) 

l yf = cf( Xi) 


where i = 1,2. Let component 1 have an adaptive key component controller with the following 
direct adaptive control law 


/ uf = Gy? 

1 G = -yf(yf) r H; H > 0 


( 8 ) 


Note that the adaptive key component controller operates only through the input and output 
ports located on component 1. The augmented control u, v (x 2 j on component 2 would be present 
if additional output feedback control were needed to make the system controllable and observable 
from the key component or to satisfy sufficient conditions for the adaptive controller, as given in 
Theorem 6.1. 

In the next section, we develop theory on the inheritance of dissipativity traits in Evolving 
Systems. This theory is used to provide sufficient conditions for adaptive key component control 
in Nonlinear Evolving Systems. 


5. Inheritance of Dissipativity Traits in Evolving Systems 

We say a trait, such as dissipativity, is naturally inherited by the Evolved System when all 
of the subsystem components have the trait and the Evolved System inherits the trait. In this 
section we show that two forms of dissipativity are naturally inherited in Evolving Systems. The 
notion of dissipativity used here is different from the dissipativity studied by other authors, e.g., 
[13, 14], 


Definition 5.1. Consider a nonlinear system of the form given by 


= A(x) + B(x)u 
- C(x) 

5 


x 

y 


( 9 ) 



We say that this system is strictly dissipative when 3V(x) > 0 Vx + 0 and V(0) = 0 such that Vx 


VVA(x) < -S (x) 

VVB(x) = C r (x) ( 

where VV = gradient V and S (x) > 0 Vx + 0. When S (x) = 0, (9) is a dissipative system. 

The function V(x) is referred to as an energy storage function or a Lyapunov function for (9). 
We can use VV to define the function V(x), as follows 

V(x) = VV[A(x) + B(x)u] (11) 

The above says that the energy storage rate is always less than the external power. This can be 
seen by using (10) to obtain 


V(x) = VV[A(x) + B(x)u] 

< —S (x) + C r (x)u (12) 

= -S (x) + (y,u) 

Taking u = 0, it is easy to see that (12) implies (10a), but not necessarily (10b). So (10) 
implies (12) but not conversely. If (10a) is an equality, then (10) and (12) are equivalent. If the 
inequalities in (10) and (12) are equalities, then the property is called strict passivity, which is 
described in Section 8. 

Definition 5.2. Consider a nonlinear system of the form given by 

(* = A(x) + B(x )u 

l y = C(x) 

We say that this system is almost strictly dissipative (ASD) when there is some output feedback, 
u = Gy + u r , that makes it strictly dissipative. 

We now give a result on the inheritance of almost strict dissipativity in systems connected in 
feedback. 


Theorem 5.3. Suppose we have a pair of subsystems of the form 


X,' 


„A _ 


where i — 1,2 and both subsystems 


Aj(Xj) + eB,(x)u, 

+ Bf (x)uf 

C,(x,) 




Cf(Xi) 




fhl 

yi 

5 

) and ( 

tel 

f y2 l 


y? 


y 2 


(14) 


are almost strictly dissipative 


(strictly dissipative) with energy storage function Vi(xi) and Vofe), and 


VV;eBj(x;) - fC' (x,) 


(15) 


Then the feedback connection of the two subsystems, where (yi = U 2 ) and (uj = — y 2 ), leaves 


the resulting composite system | u A = 
pative). 


>y a = 


almost strictly dissipative ( strictly dissi- 


I 
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Proof: We first prove the theorem for almost strict dissipativity. By the definition of almost strict 
dissipativity, there exists output feedback control that makes each of the subsystems strictly dis- 
sipative. It is well known [15], that output feedback can be added to a system (A(x), B(x), C(x)) 
to obtain 


f x = A(x) + B(x)GC(x) + B(x)u r 
1 y = C(x) 

Let uf = Gf yf + uf L where i = 1 , 2, be the output feedback that makes the subsystems given 

by (14) strictly dissipative. The subsystems with uf defined as above are now both strictly 

dissipative and can be written as 

*; = A,'(x,) + Bf (x,)Gf Cf (x,) + eB,-(x)u/ + Bf (x)uf r 

= Af(x,) + eB,(x,)u, + Bf (x ; )uf r 

y, = ckx,) <17) 

, y f = cf(x,) 

Using the definition of strict dissipativity and the assumption (15), (17) gives us 
' W/AfCx,) = Vy,[A,(x,) + Bf(x,)GfCf(x / )] 

< ~Si(Xi) 

VV/eB ,(x,) = e(C,(x,j) 7 ’ ' j 

W,eBf(x i ) - e(Cf(Xj)) r 

We can connect the two subsystems in feedback, with (y i = U 2 ) and (ui = — y 2 ). Then using 
the assumption (15), we obtain 

Vy,eB|(X|)U| = eC[(xi)[-y 2 ] 

= -eyfya 


Vy 2 eB 2 (x2)u2 = eCi, (x 2 ) [y 1 ] 


Let x = X| x 2 J , y = yf yfj and write the system from (17) in matrix form as 


A(x) + B(x)u 
Af(x x ) + eBi(xi)ui 
A, (x 2 ) + eB 2 (x 2 )u 2 
Cf( Xl ) 

C 2 (x 2 ) 


Bf (xj)l [uf 
Bf (x 2 ) uf 


Let V — Vi + V 2 . We now can use (18)-(20) to solve 

WA(x) = [TV, vv^’+^'H 
L -J |Af(x 2 ) + eB 2 (x 2 )u 2 J 

= WiAf(xi) + WieBi(xi)ui + Vy 2 Af(x 2 ) + Vy 2 eB 2 (x 2 )U2 

< -S i(xi) + (-eyfy 2 ) - S 2 (x 2 ) + ey 2 yi 

= -(S 1 (x 1 ) + S 2 (x 2 )) 

- -S(x) 
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and 


VVB(x) 


| VI/, VV 2 ] 

Cf( Xl )] r 

C,(x 2 ) 

C r (x) 


Bf( Xl ) 

0 


0 

B?(x 2 ) 


( 23 ) 


Therefore the composite system given by [ = 

0 


r 


>y a = 


output feedback 


G? 

0 


G* 


uf 

< 


is made strictly dissipative with 
Hence, the composite system is almost strictly 


dissipative. 

We now prove the theorem for strict dissipativity. Since the two subsystems are strictly 
dissipative, we can let G^ = 0 in the preceding proof. Then the proof for this result follows 
directly and the composite system is strictly dissipative. □ 

A consequence of Theorem 5.3 is that two component Nonlinear Evolving Systems with 
components that are either both almost strictly dissipative or strictly dissipative from an admittance- 
impedance point of view, inherit these properties from their subsystem components. Thus almost 
strict dissipativity and strict dissipativity are naturally inherited traits in Nonlinear Evolving Sys- 
tems. 


6. Sufficient Conditions for Adaptive Key Component Control 

We can now give sufficient conditions for the use of an adaptive key component controller to 
restore stability in Nonlinear Evolving Systems. 


Theorem 6.1. Consider a two component nonlinear time-invariant Evolving System given by 

( x, = A,(x;) + eB,(x,)u/ + Bf (x,)uf 

1 y, = c,(x,) ' 


where i = 1,2 with energy storage functions Vi(xi) and V 2 (x 2 ). Let component 1 have an 
adaptive key component controller with the following direct adaptive control law 

J< = Gy ‘ (25) 

I G = -yf(yf) r H; H > 0 (25) 

Assume that V\ and V 2 are positive Vx ^ 0, equal to 0 when x = 0 , and radially unbounded, 
that (A(x),B(x), C(x)) are continuous functions ofx, and that S (x) is positive Vx + 0 and has 
continuous partial derivatives hi x. Furthermore, assume: 

1. Component 2, given by (u 2 ,y 2 ), is strictly dissipative and in impedance form; 

2. Component 1, given by (u'^.y'j'), is almost strictly dissipative; 

3. Component 1, given by (ui,yi), is in admittance form. 

Then the adaptive key component controller given by ( 25 ) that joins component 1 with component 
2 produces global asymptotic state stability, i. e. , x = |xi x 2 J — > 0 as t — > oo with bounded 
adaptive gains and the outputs y,- = C,(x,) — > 0 as t — > oo. 
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Proof: Since the physical connection of component 1 to component 2 is equivalent to the 

feedback connection Ui = -y 2 and 112 = jq, by Theorem 5.3 we have that the closed-loop system 
(Upy|) below is almost strictly dissipative. 


*1 = Ai(xj) - eB 1 (x 1 )C 2 (x 2 ) + Bf (xi)uf 

x 2 = A 2 (x 2 ) + eB 2 (x 2 )Ci(x 1 ) (26) 

y i = cf( Xl ) 

Rewrite (25), using G[ constant to get 

uf = G iy f = G]yf + AG iy f 

AGi = G, = -yf(yf) r h 1 (27) 

where AGi = Gi — GJ and h| > 0. Combine (26) and (27) to obtain 
*1 = Af(xi) - eB 1 (x 1 )C 2 (x 2 ) + Bf (xi)wf 

x 2 = A 2 (x 2 ) + eB 2 (x 2 )Ci(xi) (28) 

yf = cf( Xl ) 

where Af(xi) = A|(X|) + Bf (xi)G,Cf (x^, u f = AGiyf, and 0 < e < 1. 

Let V — V\ + V3- We now have 


V = -5(x) + <yf,wf> 


Form Vc 


1 

2 


tr( AG 1 h j 1 AGf ) to obtain 


(29) 


V c = t^AGjh^AGf) 

= -tr(yf(yf) r h l h r i AG[) 

= -tr(yf (<yf) r ) (30) 

= — <yf ,o>f ) 

Define V(x, AG) = V(x) + Vc(AG). Then from (29) and (30), we have 

V(x,AG) s V(x) + y G (AG) 

= -S (x) + (yf , wf > - (yf , cof ) (31) 

= -S (x) < 0 


This guarantees that all trajectories (x, AG) are bounded. If V(x, AG) is uniformly continuous 
or y(x, AG) is bounded, then Barbalat’s Lemma [16] yields: S (x) — » 0, as t — > 00, and the 

positivity and continuity of S (x) imply that x = [ x 1 X 2 J — » 0 as t — » 00. 

Consider 


V(x,AG) 


= 

-S(x) 

< 

S (x) 


dS (x) . 


dx 


dS (x) 


dx 

< 

dS (x) 


dx 

dS(x) 

dx 


Pll 

[l|A(x)|| + ||B(x)|| ||wf ||] 

[l|A(x)|| + ||B(x)||||AG 1 ||||Cf(x 1 )||] 
9 


( 32 ) 


< 



which is bounded because (x, AG) is bounded, S (x) has continuous partial derivatives, and 
(A(x), B(x), C(x)) are continuous, and a continuous function of a bounded \(t) in R' v is also 
bounded in t. Therefore, y, = C,(x,) — » 0 as t — » 0. □ 


7. Special Case: Internal State Nonlinearity 


Here we look at a special case of the above theory when the only nonlinearity is in the internal 
state structure of each component. This means that the Evolving System from (7) can be written 
as 

( x; = A,x, + h,(x,) + f B,u, + Bfuf 

\ y i = C,x, (33) 

1 y t = cfx, 


where i = 1,2. The nonlinearities are assumed to satisfy h,(0) = 0 and a Lipschitz continuity 
condition given by 

||h/(x) - h f (y)|| < ju/Hx - y|| (34) 

V x, y with Hi > 0. 

We can choose quadratic Lyapunov functions of the form V (x) = \xJPjXj. We assume the 
linear part of each component (A,, B,, C,) is almost strict positive real (ASPR), i.e., 3 G* such 
that A c = A, + B,G*C, satisfies 

l 1 1 l ‘ 


/ (Af) 7 P, + P,Af = -Q,- 
\ P,B? = (C f) T 


(35) 


with P,, Q, both positive definite. In addition, we assume that 


PB, = Cj 


(36) 


with the same P, > 0. Then, from (35) - (36) we have 

( VVi-Af (x,) - xj P,-[(A/ + BfGrCfjx, + h,(x ; ) + eB,u,] 
j Q;X, + xj P,'h;(x,') + eyj U; (37) 

l Vl/,Bf(x,) - xfP,Bf - (Cfx,) 7 ' 

Now, from (34) and the Cauchy-Schwarz inequality, we have 
-lx 7 Q,x, + x 7 P/h/(X;) < 


Therefore, (37) becomes 

( VVi-Af (x,-) - x 7 P,[(A, + Bf G*Cf )x, + h,(x,) + efi,^] 

| < -S ,(Xj) + e(y„ u,) (39) 

I V\/,Bf(x,) - xfP ( Bf = (Cfxi) T 

w,th 5 ,(x,) s r,l WP > 0 V „ * 0 when r , ■ = > ft From this analysis, we 

have the following result: 


^ min m 
2 

TillXill 2 

S,(x,) 


/ ^max(^ i)H I X/ X/ 


(38) 
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Theorem 7.1. If the linear part, (A,. B.', C'.' ), of (33) is ASPR, then (36) holds, and the Lipschitz 
constant in (34) satisfies the following 


0 < Pi < 


'Pnin(Qi) 

2A max (P,) 


(40) 


Then, when component 1 is joined with component 2, the adaptive key component controller 
given by (8) maintains closed-loop stability. 

Proof: From (40), we obtain (39) and we can apply Theorem 6.1 to achieve the result. □ 


8. Inheritance of Passivity Traits in Evolving Systems 

In this section we give results describing the inheritance of strict and almost strict passivity 
in Nonlinear Evolving Systems. 

Definition 8.1. Consider a nonlinear system of the form given by 

/ x = A(x) + B(x)u 
1 y - C(x) 

We say that this system is strictly passive when 3V(x) > 0 Vx 4 0 and V(0) = 0 such that 

V(x) = (u,y> - S (x) (42) 


with S (x) > 0 Vx + 0. 


Definition 8.2. Consider a nonlinear system of the form given by 

/ x = A(x) + B(x)u 

1 y - C(x) t j 

We say that this system is almost strictly passive (ASP) when there is some output feedback, 
u = Gy + u r , that makes it strictly passive. 


We state the following result without proof. The result is an obvious consequence of Theorem 
5.3 and a well known result about passivity and systems connected in feedback [16]. 

Theorem 8.3. Suppose we have a pair of subsystems of the form 


( X, = A,(X;) + B;(X;)U; + Bf (X/)uf 

j y,- = C,(x,) (44) 

l y t = c?(x,) 


where i — 1,2 and both subsystems 


U( 


and 


u 2 


are almost strictly passive 


(strictly passive) with energy storage functions Vi(xi) and VSfe)- Then the feedback connection 
of the two subsystems, where (yi = u 2 ) and (ui = — y 2 ), leaves the resulting composite system 


u A = 


■ y a = 


almost strictly passive (strictly passive). 


11 



Theorem 8.3 can be used to show that two component Nonlinear Evolving Systems with 
components that are either both almost strictly passive or strictly passive from an admittance- 
impedance point of view, inherit these properties from their subsystem components. Thus strict 
passivity and almost strict passivity are naturally inherited traits in Nonlinear Evolving Systems. 

All of the preceding results assume that the Lyapunov function, V(x), is defined on the entire 
domain, R”, of the system. Thus all the stability and dissipativity results and Theorem 6.1 are 
global results. However, the Lyapunov function, V(x), might only be defined on a neighborhood 

IV, (0, rf) s i — < r j l of the origin, in which case the results could only be local at the best. 

lllx/ll J 

Using Lemma 1 from [17], 35 > 0 such that if the initial conditions of the system are close 
enough to the origin, i.e., within Ng = (0,5), then the trajectories are guaranteed to stay in the 
neighborhood of the origin for which the Lyapunov function is defined. In such a case, then the 
results would be local. Lor instance, if the Lyapunov function L(x) in Theorem 6.1 only has the 

assumed properties on a neighborhood IV, (0 r,T of the origin and the trajectories all 

remain inside the neighborhood, then the stability is locally asymptotic to the origin. In that case. 
Theorem 6.1 gives the result that a Nonlinear Evolving System with an adaptive key component 
controller as given by (25) will have bounded gains and locally asymptotic state tracking. 


9. Equivalence of Traits in LTI Evolving Systems 

In this section we give equivalence results for some properties of LTI Evolving Systems. 

Theorem 9.1. If a nonlinear system given by (A(x), B(x), C(x)) is strictly passive, then it is 
strictly dissipative. 

Proof: Let V(x) be the energy storage function for the system (A(x), B(x), C(x)). By the defini- 
tion of the gradient and strict passivity, we have 

V(x ) = VVA(x) + VB(x)u 
= -S (x) + y r u 


Let uiO. Then (45) becomes 


V(x) = VVA(x) = -S (x) < -S (x) (46) 

Thus, y r u = WB(x)u Vu. Since y = C(x), this implies that WB(x) = C r (x). Therefore, the 
system (A(x), B(x), C(x)) is strictly dissipative. □ 

Theorem 9.2. A LTI system given by (A, B, C) is strictly dissipative iff it is strictly passive. 

Proof: Lirst we show that a LTI system that is strictly dissipative is strictly passive. Since 

(A, B, C) is LTI and strictly dissipative, by the Kalman-Yacubovic Lemma [16], 3 quadratic 
Lyapunov function V (x) such that V x e R“ 

VVAx = ix r (PA + A r P )x<-S(x) 

VLB = x r PB = x r C = y r (4/) 
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Choose S (x) = ^x r Qx with Q > 0. Then we have 

VVAx = ±x r (PA + A r P)x < -±x r Qx 

= PA + A r P < — Q “ (48) 

Let W = PA + A 7 P, then W is symmetric and W < -Q. Then 3H, H symmetric such that 
W + H = -Q. This can be seen if we let H = -(Q + W) > 0. Then Q + H > 0, because Q > 0. 
So 

VVAx = PA + A r P = -(Q + H) (49) 

Now we can write 

V(x) = WAx + WBu = -(Q + H) + y r u (50) 


V x e R". Hence, (A, B, C) is strictly passive. The converse is true by Theorem 9.1. □ 

Theorem 9.3. A LTI system given by (A, B, C) is strict positive real ( SPR ) iff it is strictly passive. 


Proof: First we show that if (A, B, C) is SPR, then it is strictly passive. Since (A, B, C) is SPR, 
the Kalman- Yacubovic Lemma [16] implies that 3e>0 such that 


( (A + eI) r P + P(A + el) = -Q 
\ PB = C T 


with Q > 0 and P > 0. We can rearrange (51) to obtain 

I A r P + PA = -(Q + 2eP) 
\ PB = C T 


(52) 


Since P > 0 and Q > 0, then W(e) eQ + 2eP > 0. Choose V(x) s ^x r Px with P chosen as in 
(51) and P > 0. The time derivative along any state trajectory of L(x) is given by 

V(x) = \ (x 7 PAx + x r PBux) 

= |x r (A r P + PA)x + x r PBu 
= -jX r W(e)x + (Cx) r u 
= -S (x) + y r u 

Therefore (A, B, C) is strictly passive. 

We now show that if (A, B, C) is strictly passive, then it is SPR. Since (A, B, C) is strictly 
passive, we have that V(x) = -S (x) + y 7 ti with S (x) > 0. Choose S (x) = W(e), with P and Q as 
defined in (51). Then all of the previous arguments can be reversed, giving the desired result. 

We have shown that for LTI systems, strict positive real is equivalent to strict passivity. □ 

A system that is almost strict positive real (ASPR) is defined to be one that can be made strict 
positive real with output feedback. Hence, for LTI systems, the above results give us that almost 
strict positive real is equivalent to almost strict passivity. Likewise, the above results give us that 
almost strict passivity is equivalent to almost strict dissipativity in LTI systems. 

The equivalence results of the above properties in LTI systems also apply to LTI Evolving 
Systems. Hence, for LTI Evolving Systems, we can substitute any of the equivalent properties 
into Theorem 6.1. 
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10. Illustrative Example of Nonlinear Evolving System with an Adaptive Key Component 
Controller 


We now give an illustrative example of an adaptive key component controller restoring sta- 
bility in a Nonlinear Evolving System where each of the components is actively controlled and 
stable, but the Evolving System fails to inherit the stability of the components. Consider a fully 
actuated, fully sensed three mass Evolving System. Component 1 has only one mass. The dy- 
namical equations for component 1 are 


1 rnqi 

= U\ 

1 * 

= [<7i <7i] 


(54) 


where m \ = 30 is the mass of mass 1, q\ is the displacement of mass 1, and u\ = -(0.9s + 0.1)gi 
is the local controller for component 1, where s is the Laplace variable. 

Component 2 has two masses connected by a nonlinear spring. The dynamical equations for 
component 2 are 


m 2 q 2 = u 2 - k 2i (q 2 - qi) 

mh = M2 - [&23(<?3 - qi) + yU23 sin(g 3 - q 2 )] 

yi = [<?2 < 72 ] 

. ys = [<73 <73] 


(55) 


where m 2 = 1.0 is the mass of mass 2, m 2 = 1.0 is the mass of mass 3, q 2 is the displacement of 
mass 2, q 2 is the displacement of mass 3, k 22> - 1.0, and n 22 — 0.5. The controllers on component 
2 are 


/0.1 \ 
u 2 — — I — + 0.2s -(- 0.5 1 q 2 

m 3 = - (0.6s + 1) q 2 


(56) 


The controllers for components 1 and 2 have been designed to produce stable behavior when the 
components are unconnected. The two components are joined by a spring, k\ 2 = 1.0 connecting 
mass 1 with mass 2. The connection parameter, e, multiplies the spring connecting the two 
components. 

A Simulink model of this system has been created to study its stability. The Evolving System 
was evaluated in its fully evolved state, where the local component controllers were connected 
to their inputs and outputs and e was set equal to 1 . Figure 1 shows the nondimensional mass 
displacements for the Evolving System given by (54)-(56). When the system is fully evolved, 
i.e., e = 1, the Evolved System is unstable as seen in Fig. 1. 

A Simulink model was created to implement the adaptive key component controller given 
by (8) for the illustrative Nonlinear Evolving System example. Simulations were run with the 
connection parameter, e, ranging from 0 to 1, allowing analysis of the full evolution of the system. 
The key component controller was able to maintain system stability during the entire evolution 
process when it used the input-output ports on mass 1 of component 1, see Fig. 2. 


11. Conclusions 

We provided the reader with the framework for Nonlinear Evolving Systems of flexible struc- 
tures. We developed the idea of modeling the contact dynamics between components as the feed- 
back connection of the impedance of one component with the admittance of a second component. 
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A major contribution made by this paper was the use of the impedance-admittance formulation to 
deliver theory on the inheritance of dissipativity, passivity, and positive real traits of systems con- 
nected in feedback. Using the impedance-admittance formulation of contact dynamics between 
components of Evolving Systems, we were able to show that strict and almost strict dissipativ- 
ity and passivity traits are naturally inherited in Nonlinear Evolving Systems. We also showed 
that for LTI systems, strict dissipativity, strict passivity, and strict positive realness are equivalent 
properties, with the strict dissipativity trait having the weakest hypotheses. A principal contri- 
bution of this paper is the development of the adaptive key component controller and sufficient 
conditions for its use in Linear and Nonlinear Evolving Systems. Finally, we demonstrated the 
adaptive key component controller restoring stability in an illustrative example of a Nonlinear 
Evolving System. 
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Figure 1: Nondimen sional mass displacements of example with no key component control. 



Time (seconds) 


Figure 2: Nondimensional mass displacements of example with adaptive key component controller on component 1. 
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